RADIX AND PSEUDODIGIT REPRESENTATIONS IN Z 



EVA CURRY 

Abstract. We define radix representations for vectors in Z™ analogously with radix representa- 
tions in Z, and give a sufficient condition for a matrix A : Z n — ¥ Z n to yield a radix representation 
with a given canonical digit set. We relate our results to a sufficient condition given recently by 
Jeong [2] . We also show that any expanding matrix A : Z" — ¥ 1 n will not be too far from yielding 
a radix representation, in that we can partition J7 1 into a finite number of sets such that A yields 
a radix representation on each set up to translation by A N s for some vector s (N > will vary). 
We call the vectors s pseudodigits, and call this decomposition of Z™ a pseudodigit representation. 



1. Introduction 

A number x in some set (such as the whole numbers, integers, or real numbers) is written in a 
radix representation if there are a base (also called a radix) j3, a finite set of digits D, and a finite 
number N such that x can be written 

N 

3=0 

with the dj e D. Examples of radix representations for the whole numbers include our standard 
decimal representation as well as binary representation. In both these examples, the number is 
usually written as a string of digits djvdjv-i* • • dido, where the position of a digit in the string 
indicates the power of j3 that the digit is multiplied by. Thus radix representations are also known 
as positional number systems. It is well known that the nonnegative integers can be written in radix 
representation with any positive integer /3 > 2 as a base and with digit set [/?] = {0, 1, . . . , j3 — 1}. 
Matula [6] has shown that digits sets can be found to write the integers (both positive and negative) 
in radix representation if the base satisfies \/3\ > 3, (3 6 Z. 

Radix representations, or positional number systems, can be generalized to higher dimensions. 
For arbitrary dimension n > 1, any matrix with integer entries A £ M„(Z) defines a linear trans- 
formation on Z n . We consider expanding matrices in the following sense. 

Definition 1.1. A dilation matrix is a matrix A with integer entries, each of whose eigenvalues A 
satisfies |A| > 1. 

The determinant of a dilation matrix is an integer q = |detA| with q > 2. A dilation matrix 
gives a mapping of Z™ into a sublattice of Z™ with nontrivial cokernelQ 
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Let A be a dilation matrix and let D be a set of coset representatives for Z"/yl(Z"). Note that 
D has q elements. Any x G Z" is in a unique coset of Z ra /A(Z n ), and thus can be written uniquely 
as 

x = Ay + r, 

for some y £ Z™ and r E D. This defines a Euclidean algorithm for decomposing any xo S Z": 

.t = An + r 
xi = Ar 2 + ri 
a?2 = ^£3 + r 2 

where the Xj and r$ are unique, given xq. We say that this algorithm terminates if there exists 
some N > 1 such that xj = rj = for all j > N. 

Definition 1.2. An n x n matrix A yields a radix representation (for Z") with digit set D if the 
Euclidean algorithm described above terminates for every x G Z™. 

That is, if A yields a radix representation with digit set D, then for every x £ Z n there exists a 
nonncgative integer A = N(x) such that 

AT 

x = Aiei + d = ^ 2 ^2 + ^4di + d = ■ ■ ■ = A j d 

3=0 

for some sequence of digits d 3 ; e D. Thus we can represent x by the sequence 

x = (d n , . . .,dx,d ), 

just like positional number systems in one dimension. As in the case of positional number systems 
for Z, we call the coset representatives dj digits, and to call D a digit set. 

Not all dilation matrices yield radix representations. For example, when A = 2 (n = 1), with 
digit set D = {0, 1} we can represent every positive integer with a radix representation (writing 
them in binary), but we cannot represent any negative integers. When A docs not yield a radix 
representation, we can still find a unique sequence of digits associated to each x 6 Z™ using the 
algorithm above. However this sequence will be infinite for vectors for which the algorithm does not 
terminate. Section [5] presents an alternative representation in the case that a dilation matrix does 
not yield a radix representation. In section [3l we give sufficient conditions for a dilation matrix A 
to yield a radix representation. 

The existence of a radix representation depends on the choice of digit set as well as the matrix 
used as a radix. A digit set can be any set of coset representatives, so many potential digit sets 
are possible. In general, if one takes any fundamental domain F for Z" (that is, any set F C K ra 
congruent to M"/Z"), then A(F) n Z" will be a set of coset representatives of Z"/^l(Z n ), and thus 
a candidate for a digit set. For positional notation in one dimension, the digit set is usually taken 
to be {0,1,...,|/3|}. This is not a satisfactory digit set in the present context, since it will not 
allow us to represent negative numbers with a radix representation when /3 is positive. We run into 
problems since the origin is on the "edge" of the digit set in some sense (we will make this more 
precise in section |3]). In general, choosing the smallest possible coset representatives as digits, and 
choosing digits that can be generated by the intersection of a convex set with Z 71 , seems to produce 
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digit sets that are easier to work with and more likely to yield a radix representation. Thus in this 
paper we consider the digit set 



While not all dilation matrices give a radix representation, we give a representation that always 
holds. 

Definition 2.1. When an n x n dilation matrix A does not yield a radix representation with 
digit set D = A(F) n Z n , but does satisfy Theorem 12.21 below, we say that A yields a pseudodigit 
representation. 

Theorem 2.2. Let A be a dilation matrix of dimension n > 1. Let D — A(F) l~l Z n , where 
F = [— 2> o)™' There exists a finite set S such that every vector x 6 Z™ can be written in exactly 
one of the following forms: 

(1) x — ^2j=o A^dj , where dj G D for each j and N(x) > 0; 

(2) x = A N( > x h + ES'" 1 AH], where d 3 G D, s G S, and N(x) > 0. 

That is, Z 71 can be partitioned into a set of vectors that can be written in a radix representation, 
and (a finite number of) translates of that set by sets {A^s : j > 0, s G S}. Call the elements of S 
pseudodigits. 

Let us clarify the notation. Individual vectors x G Z™ may have a radix representation with 
matrix A and digit set D, however we only say that A yields a radix representation if every x G Z™ 
has a radix representation (and thus no x has only a pseudodigit representation). Those vectors 
which do not have a radix representation will always have a pseudodigit representation instead. To 
distinguish between the case where all vectors x have a radix representation, and the case where 
some vectors x have only a pseudodigit representation instead of a radix representation, we say 
that A yields a radix representation with digit set D, or A yields a pseudodigit representation, 
respectively. We will also say that a vector ieZ™ has a pseudodigit representation only when x 
does not have a radix representation. 

We must first define a norm under which the dilation matrix A is expanding. The I 2 norm is not 
sufficient for our purposes, since ||Ae||2 may not be greater than ||a;||2 for all vectors x G M™. Take, 
for example, the dilation matrix and vector 



D := A(F)f]Z n , 



where F = [-§, ±) 



2. Pseudodigit Representation 



A = 



2 -2 
1 2 



x = 



1 
1 



Lemma 2.3. Let || • || 



' be the norm defined by 




Then \\Ax\\' > \\x\\' for all f x G Z n . 



Proof. We can see that this series converges using the root test and the fact that the spectral radius 

r(B) = max{|A| : A an eigenvalue of B} 
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for any matrix B can also be found by the formula [5] 

r{B) = lim 

j-s-oo 

where ||-B|| op is the operator norm 

||B|| op = sup{||Ba;|| 2 : x G R", ||x|| 2 = 1}. 

Since 

(\\A- j x\\j) l /i < (WA-iWD^iWxg) 1 ", 
with limj-voodla;!!!) 1 ^ = 1 for all x G R™, 

lim (IIA-^HDVj < max{|A|- 2 : A an eigenvalue of A} < 1, 

so our series does converge, and defines a norm on R™. Note that, for x ^ 0, 

/ oo \ 1/2 / oo N 1/2 

\\Ax\\' = lj2\\A- j+1 x\\j) = ll^" 2 

X 1/2 

|Ac||2 + £||A-^||l >|N|'. 

□ 

Remark 2.4. Our norm is similar to the || • ||' norm used by Lagarias and Wang [3]. However, they 
use 

oo 
3=0 

where p may be any scalar satisfying 1 < p < min |A|, with the minimum taken over all eigenvalues 
A of A. We use the I 2 version of this norm, with p = 1, since we will often wish to compare ||x||' 
and ||x||2 for vectors x € 1". In particular, while the traditional definition of a dilation matrix 
involves a lower bound on the modulus of the eigenvalues, we will be considering lower bounds on 
the singular values, in which case it will be more natural to consider the I 2 norm. 



Proof of TheoremUTM Let 

m : = min {|| Ar||' : ||ir||' = 1}, 

M : = max{||Ar||' : ||x||' = 1}, and 

p:=max{||/||': / e F}. 

Then ||r||' < Mp for every digit r E D. 

Given x £ Z n , let y e Z" and r G D be the unique vectors such that x = Ay + r. Then by the 
definitions of m, M , and p, 

I,.,,/ ii,-i/ mi// ii / NI' + IHI' ^ \\x\\' + m p 

y =\\A {x-r)\\ <— \\x-r\\ < < 

m mm 

If x is big enough, \\y\\' will be less than ||x||'. A sufficient condition for this is 

M> \\^l±MP i 
m 



RADIX AND PSEUDODIGIT REPRESENTATIONS IN Z' 



5 



which simplifies to 

Mp 



> 



m — 1 

That is, the sequence (xq, x\, X2, • • •) generated by the multidimensional Euclidean algorithm de- 
scribed above is decreasing until some xn is in the ball of radius R :— -^3[- Once inside this ball, 
the subsequent Xj (for j > N) remain in the ball, since in each case 1 1 ajj 1 1 1 ' < — jr and thus 

+ Mp M p 



< 



1 



For some vectors xq, the sequence (xj) will continue to decrease down to the origin. That is, we 
can find a (finite) radix representation for these vectors. For vectors xq that do not have a radix 
representation, since the Xj take values in a finite set for j > N, there will be some k > 1 and some 
I > N such that Xj+k = xj for all j > I (the sequence will start to repeat at some stage after it 
has entered the ball of radius R). There are a finite number of possible cycles that the sequence 
(xj) can fall into, again since there are a finite number of integer vectors inside the ball of radius 
R. Choose one vector in the cycle to represent each distinct cycle. These cycle representatives are 
the pseudodigits. □ 

This suggests a method for determining whether a given dilation matrix A gives a radix represen- 
tation for every ieZ™. When A is normal, we can write x = x\, where the sum is taken over 
all eigenvalues A of A, and x\ is the orthogonal projection of x onto the eigenspace corresponding 
to the eigenvalue A. We may then write 

oo oo oo 

Nir 2 = X>*H' = £E \\A- j x x \\j = J2J2 u-'xxwi = £ EaApnwi* 

A A j=0 A j=0 A j=0 



and 



A 1 ' A 

'2 ^ i\/t2\ 



< M 2 1 1 x 1 1' 2 , with (since A is normal), 

m = min{|A| : |A| an eigenvalue of A}, 

M = max{|A| : |A| an eigenvalue of ^4}. 

We can then compute the radius R and check whether A gives a radix representation of each vector 
x e Z™ with ||x||' < R. 

If the dilation matrix A is not normal, it may be difficult to find R. We may use the singular 
values 8 of A to compute an upper bound for the radius i?, however, provided a > 1 for every 
singular value a of A. Let 

/i = min {er : er a singular value of A}, 
v = max {er : a a singular value of ^4}. 

If p > 1 then 

||Ar||' 2 > IIArll 2 : > /i 2 !!*!! 2 . 
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for all x £ M. n , so m > \i. Note also that fi 1 is the maximum of the singular values of A 1 , with 
fi- 1 < 1. Then || J 4"-'a;||| < Gti 2 ) _J '||a;||| for all x G K", and 



\\Ax\f < ^)- j \\M\l = J-rWMl < -t-t- 2 M 2 

3=0 

(2 2 \ I/ 2 
-^j • Now we may give a bound on R 

m 



.■/2 

v/n 



i? < 



2(/i-l) 2( M - 1)3/2^ + 1)1/2- 

A program for finding digits and pseudodigits for dilation matrices can be obtained by contacting 
the author. A few specific examples are given in sections [4] 

3. Radix Representation 

As above, set 

fi = min{(7 : a a singular value of A}, 
v = max {a : a a singular value of A}. 

If /i is large enough, in particular, if /i > 2y / n, then n Z™ C £>, where is the ball of radius 
•y/n in the || • ||2 norm centered at the origin. This means that all vectors of the form 



o 2 



4,0,1 for i = 1,2, 



are contained in the digit set D. As the following theorem shows, this gives a sufficient mathematical 
formulation of the intuitive idea that we need the origin to not be on the "edge" of the digit set in 
order to ensure that A yields a radix representation. 

Theorem 3.1. Let A be an n x n dilation matrix. If 

[x > 2\fn 

then A yields a radix representation of Z™ with the canonical digit set D = AF n Z™, where F = 

I 2' 2/ 

Proof. Let xq be a vector in Z n and let {xj : j > 0} be the orbit of xq under repeated applications 
of the Euclidean algorithm. That is, xo = Ax\ + tq, x\ = Axi + r*i, and so on, with the rj in our 
digit set for j > 0. We showed in Theorem 12.21 that there exists some N = N(xo) > 1 such that 
either Xj = r j = for all j > N, or xn is a pseudodigit. Recall that when xn is a pseudodigit, the 
orbit of Xo under the Euclidean algorithm falls into a cycle that includes xn, so that there exists 
an integer I > 1 and digits fo, • . • , rj_ i such that 

l— l 

i=0 

We use our canonical digit set D — A(F) n Z ra , so that for each j > 0, rj = d for some d E D. 
Given any a; € Z™ and > 1 as above, if Xj — for every j > N, then we have found a radix 
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representation for x . Thus suppose that £ at is a pseudodigit. Write x in place of a; at, since the 
subscript is no longer necessary. Note that x n = A N ~ 1 x + d for some d € Da,n, where 

3 

D AJ = {£ A l di : di e D} for any j > 1. 

i=0 

Thus xq has a radix representation if and only if x has a radix representation. 

When x is a pseudodigit, there exists an integer I > 1 and a vector d x G Da,« such that 

We rewrite this as 

A~ e x = x + A~ e d x 
x = A~ e x + g x 

for g x = —A~ l d x . Then 

x = g x + A~ e x 
= g x + A- f g x + A~ n x 



= J2A-Vg x =J2(At)-ig x . 
3=1 j=i 

We can now bound the size of x as follows (since A~ l; >x — > as j ' ^ oo). 

00 00 / 1 V 1 

lkl| 2 < E IK^)"^x|| 2 < E 77 ll^lla = rll^lh, 

3 = 1 3 = 1 VM 7 M 



with 



£-1 



|.9,|| 2 - II - .g*|| 2 < £ ||A-^-|| 2 = £ Il^- J "^ll2 

j=0 j' = l 

<Ell^.M| 2 for some / j6 F 



i=0 
i=0 



1 yri [i—l [i \fn 



1 [1—1 [1 y/n 



since D = A(F) n Z n . Thus 

" X " 2 -^-l [i e [i-12 ~ 2[i^{[i-iy 
Note that 

1 1 

< 



2[i ( - 1 ([i- 1) " 2 

when jU > 2y / n > 2, thus a; € -By^nZ™ C D. That is, a; is in fact a digit and not a pseudodigit, and 
xq has a radix representation. Since our choice of xo was arbitrary, A yields a radix representation 
for every x G Z™, with digit set D. □ 
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Remark 3.2. A similar calculation shows that ||x||' < y/n as well; however the I 2 bounds are 
enough to show our result. We could also replace the hypothesis fi > 2y/n with the slightly weaker 
hypothesis to > 2 v / n and follow through the same steps as above using the || • ||' norm instead of 
the || • || 2 norm. We chose to state the theorem with the hypothesis on /i since the singular values 
of a matrix can easily be computed, whereas checking if to > 2-y/n may be computationally more 
difficult. 

We note that any dilation matrix A is eventually expanding in the I 2 norm, in the following 
sense. 

Corollary 3.3. For every dilation matrix A, there exists a positive integer (3 > 1 such that A@ 
yields a radix representation with digit set Dp = A@(F) n Z". 

Proof. In the proof of Theorem 12.21 we saw that 

oo 
3=0 

converges (using the root test and the spectral radius formula). This implies that 

lim ||A- J || op = 0. 

However, 

\\A~ J \\op = max {cr : a a singular value of A~ 3 } 
= (min{cr : a a singular value of A° \) 

for each j > 0. Thus, setting 

fj,j = min{<7 : a & singular value of A J }, 

for each j > 1 , we see that 

lim /ij = oo. 

j-s-oo 

In particular, there exists an integer j3 > 1 such that 

Hp > 2\fn. 

We then apply Theorem 13. II to the dilation matrix A@ . □ 

By the corollary, we see that we can recapture and improve upon Theorem 13.11 using a result of 
Jeong *2 that the author has recently become aware of. Jeong defines Co to be the set of vectors 

ei 

(-2 

such that ei = for all but one index i — io, with Ci = ±1. We then set 

C = C Q {J{0} = B 1 f)Z n . 

Theorem 3.4 (Jeong). Let M be an n x n matrix with integer entries and U := (— i, 5]™- If 

(1) C C MU (or equivalently C C MU D Z n ), and 

(2) lim^ocM^ = W\ 
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then M yields a radix representation of Z™ . 
Corollary 3.5. Let A be an n x n dilation matrix. If 

(i > 2 

then A yields a radix representation of Z™ with the canonical digit set D = AFC) Z™, where F = 

L 2' 2^ ' 

Proof. Note that using Jeong's fundamental domain [/ is equivalent to using our fundamental 
domain F. If A is a dilation matrix, then lim^oo A'F = M™ by Corollary 13.31 If fi > 2, then 
Si C ABi/ 2 . Since £?i/ 2 C F (equivalently, B 1/2 C U) and C C Si, this implies that C C 
(equivalently, C C At/). The result then follows from Theorem 13.41 □ 

Remark 3.6. Consider the dilation matrix A = diag{2, 2, . . . , 2}, the nx n matrix with every 
diagonal element 2 and every non-diagonal element 0. As in the first example in the following 
section, this matrix does not yield a radix representation of Z™. However, for this matrix /i = 2, 
and indeed every eigenvalue A = 2. Thus the sufficient condition fi > 2 in Corollarv l3.5l is sharp. 

Remark 3.7. Unlike in the proof of Theorem 13. 1[ we cannot use the || • ||' norm to replace the 
condition fi > 2 with the same condition on the eigenvalues of the dilation matrix A. 



4. Examples 

4.1. Pseudodigit Representations in One Dimension. When n — 1 and A — 2, the digit set 
D = A(F) n Z is 

D= [-1,1) nz = {0,-1}. 

Note that this is not the standard digit set D — {0, 1} for binary numbers, so we are able to 
represent negative integers rather than positive integers with a radix representation. In this case, 
m = M = 2 and R = 1. Thus we only need check the set {—1,0,1} to find the set S of pseudodigits 
for A. It is easy to see that —1 and have radix representations. 

-1 = 2°(-l) 
= 2°(0) 

However, if we apply the Euclidean algorithm to 1, we find that 

l = 2(l) + (-l) 
l = 2(l) + (-l) 



The algorithm does not terminate, but repeats through the cycle (1) (of length one). Thus we 
have one pseudodigit. Since 1 is the only element in this cycle, we will take 1 to be our cycle 
representative. Then S = {!}, and A = 2 yields a pseudodigit representation. 
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Figure 1. A{F) in two dimensions 



4.2. Radix Representations in One Dimension. If we take A = —2, then again D = {0, — 1} 
and R — 1. In this case, however, every x within the ball of radius R = 1 centered at has a radix 
representation: 

-l = (-2)°(-l) 

= (-2)°(0) 

1 = (-2) 1 (-l) + (-2)°(-l) 



Thus A = — 2 yields a radix representation. 

When A — /3 <E Z with |/3| > 3, Theorem 12.21 shows that A yields a radix representation for Z 
with our standard digit set D. 



4.3. The Twin Dragon. The twin dragon matrix is the matris0 

A = 



1 1 

-1 1 



It is one of the smallest examples of an expanding 2x2 matrix, and dilation by the twin dragon is 
analogous to dilation by 2 in many ways. The set A([— \, |) 2 ) is shown in figured] 
The digit set D for the twin dragon is the intersection of this set with Z 2 : 



-1 




}• 



The twin dragon matrix gets its name from the shape of the set 



T = {V\l ■'•/, : G D}. 
3=1 



For more information, see [5]. 
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(Note that det A = 2 = \D\.) As in the case of dilation by 2, we have only two digits, and we need 
one pseudodigit for a pseudodigit representation of all vectors in Z 2 : 


u 1 

We can then partition Z 2 into "negative" and "positive" vectors (those with radix representations 
and those with pseudodigit representations, respectively), as in the case of dilation by 2. Figure [5] 
shows the set of points representable with a radix representation J2j=o -^dj with N < 6 and dj E D 



together with the set of points representable with a pseudodigit representation A n s 



with N < 6 and Dj e D. 

In order to satisfy the sufficient condition for a radix representation of Z 2 , we must square the 
twin dragon matrix: 

2 



B = A Z = 



-2 



Then detB = 4 and the eigenvalues of B are Ai = iy/2 and A2 = —iV2, with | A x | = |Aa| = V2- 



4.4. A Higher Dimensional Example. Another dilation matrix of interest in the study of mul- 
tivariable wavelets is the following matrix]^ 



A = 






1 














1 











-1 


2 


1 





-1 


1 



^This matrix, found by Lagarias and Wang [4], is the smallest known example of a dilation matrix for which a 
Haar-like scaling function for a wavelet set cannot be constructed. 



-10 



z ■ = ■ i 

■ i i 

:!!■ 
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lull 
ii 1 1 

i j i i = 



l,il: 



Figure 2. Radix and pseudodigit representations with N = 6. Boxed point are 
those with a radix representation. Diamond points are those with a pseudodigit 
representation. 
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This matrix has determinant detA = 2 and eigenvalues Ai = — yJ—\ — A2 = \ + 

h + and A 4 = \ — \ — The digit set for this matrix is 



" " 




" 





,di = 








-1 


_ _ 




-1 


" " 




" -1 " 


1 










,92 = 














A3 = \/ —f + "2" ' aI1U Aa ~ V 2 



D = {d = 

In this example, there are two pseudodigits 

S = { Sl = 



5. General Lattices 

Any n-dimensional point lattice T in K" can be expressed as M(Z n ) for some matrix M € M n (R), 
where M has full rank, and is thus invertible pQ. That is, any n-dimensional lattice T in K™ is 
isomorphic to Z™. Similarly, we can consider point lattices T in C", where T = LIT 1 for some 
invertible matrix L £ M n (C), and thus is also isomorphic to Z™. 

Let T be an n-dimensional point lattice in IR ra or in C n , as above. Let e\, . . . , e n be the canonical 
basis for Z": for each i = l,...,n, e, has I*' 1 entry 1 and all other entries 0. Then fi := Mei (or 
/i := Lei), i = 1, ... ,n, is a basis for T. For each a; in T, we may write 

a: = Xifi + X2/2 H h ain/n. 

for some coefficients 1, eZ, ! = l,...,n. Under the norm 

/ n \ V2 
|N| i3(r) := (EM 2 ) ' 

the lattice T is isometrically isomorphic to Z n . 

We call an n x n matrix A a dilation matrix for the lattice T if it gives a mapping A : F — > F 
with nontrivial cokernel. That is, A is a dilation matrix for r if A = MB for some dilation matrix 
B (for Z"), where M is the matrix giving the isomorphism between Z™ and F, V = Mil 1 . Notice 
that A yields a radix representation of T if and only if B yields a radix representation of Z™ . The 
digit set for this radix representation can be written as T H A(F) for some fundamental domain F 
for the lattice F. Define Fp to be the specific fundamental domain 

Fv = MF = M{[- 1 -, 

Let A be a dilation matrix for the lattice T. Corresponding to the norm 
define the norm 



j' on Z™, we may 



00 

E 

3=0 



A~ 



U 2 (r) 



on r. We can now restate our results in greater generality. 

Theorem 5.1 (Pseudodigit Representation). Let T — MZ n be a lattice in R™ or C™. Let A be a 
dilation matrix for F. Let D — A(Fy) H T. There exists a finite set S C T such that every vector 
x G r can be written in exactly one of the following forms: 
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(1) x = J2j=o A^dj, where dj e D for each j and N(x) > 0; 

(2) x = A N ( x h + E2o )_1 Ai *i> where d 3 ^D, seS, and N(x) > 0. 
Let 

// := min{||Aa;||;2(r) : x € T, ||a;||;2(r) = 1}. 
Note that /x is equal to the smallest singular value of B — M~ 1 A (where B is a dilation matrix for 
Z"). 

Theorem 5.2 (Radix Representation). Let T = MZ™ be a lattice in R™ or C". Let A be a dilation 
matrix for V. If fi' > 2 then A yields a radix representation of T with digit set D = A(Fr) H T. 

As a result of the Radix Representation Theorem, we also have the following result. 

Lemma 5.3. For every dilation matrix for T, A, there exists a positive integer /3 > 1 such that A@ 
yields a radix representation ofT with digit set Dp = A l3 (Fr) H T. 
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